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Abstract. The problem of reconstruction a function from spherical means is at the heart of 
several modern imaging modalities and other applications. In this paper we derive universal back- 
projection type reconstruction formulas for recovering a function in arbitrary dimension from averages 
over spheres centered on the boundary an arbitrarily shaped smooth convex domain. Provided that 
the unknown function is supported inside that domain, the derived formulas recover the unknown 
function up to an explicitly computed smoothing integral operator. For elliptical domains the integral 
operator is shown to vanish and hence we establish exact inversion formulas for recovering a function 
from spherical means centered on the boundary of elliptical domains in arbitrary dimension. 
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1. Introduction. Let SI C M. n be a bounded convex domain in M. n with smooth 
boundary. In this paper we study the problem of recovering a function / : R" — > M. 
with support inside fl from the averages (spherical means) 



over spherical surfaces with centers x £ dil and radii r > 0. Here S n ~ 1 C R ra is the 
n— 1 dimensional unit sphere, w„_i its total surface area and dS is the standard surface 
measure. We develop closed form reconstruction formulas of the backprojection type 
for recovering the function / from its spherical means Aif (x, r) defined by The 
derived formulas can be applied to arbitrarily shaped domains in arbitrary dimensions 
and recover the unknown function modulo an explicitly computed smoothing integral 
operator K,q. For elliptical domains, the operator ICq is shown to vanish. We therefore 
establish exact reconstruction formulas of the backprojection type in these cases. Our 
results generalize the ones recently obtained in [10] for n = 2 and in [14] for n = 3 to 
the case of arbitrary spatial dimension. 

The problem of recovering a function from spherical means is at the heart of 
many modern imaging applications, where the centers of the spheres of integration 
correspond to admissible locations of detectors recording some physical quantity en- 
coded in /; see Figure ITTTI For example, recovering a function from spherical means 
is essential for the hybrid imaging techniques photoacoustic tomography (PAT) and 
thermoacoustic tomography (TAT) where the function / models the initial pressure of 
the acoustic field induced by a short electromagnetic pulse. In these applications the 
inversion from spherical means arises in three spatial dimensions (see [9l 111! [24] ) as 
well as in two spatial dimensions in variants of PAT/TAT using integrating detectors 
(see [3 HU [25]) instead of the more common point like detectors. In fact these ap- 
plications initiated the authors interest in the problem of recovering a function from 
spherical means. The inversion from spherical means is, however, is also essential for 
other technologies such as SONAR (see [2j[21]), SAR imaging (see [H[22]), ultrasound 
tomography (see [H1[T7]), or seismic imaging (see [3] [6]). 
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{xi G R" : \x-xt\ = r} 



Fig. 1.1. Recovering a function from spherical means. Suppose that the function f (representing 
some physical quantity of interest) is supported inside the domain f2. Detectors are placed at various 
locations x on the boundary dQ of the domain and record averages of f over spherical surfaces 
{xi S K n : \x — x\\ = r} with radii r > 0. In this paper we derive explicit formulas for recovering 
the function f from these spherical averages (see Theorems \1.2\ [731 and \l.4\ l- 



1.1. Main results. Before we present our main results, we introduce some no- 
tation. For any integrable function ip : M™ — > R, we define the Radon transform 



(TZip) (w,s) 



<p (suj + y) dS (y) for (w,s)gS'™" 1 x 



where uj 1 - :— {y G 1" : u ■ y — 0} denotes the hyperplane consisting of all vectors 
orthogonal to lu G S n ~ 1 . The derivative of a function ip : S 1 ™ -1 x R — )• R in the second 
argument will be denoted by (d s ijj) (w, s), and (H s ip) (u, s) will be used to denote 
the Hilbert transform in the second argument (defined as the convolution with the 
principal value distribution 1/(t:s)). Further, for two distinct points xq,x± G R™, we 
set 



to* (x ,xi) := 



x\ - x 
\xi - x \ 



s+ (x ,xx) 



1 \xi\ 2 - \x \ 2 

2 |a;i — ac 1 



(1.2) 



As illustrated in Figure fOl the set H* (xq, xi) = {x G R" : oj± (x 0l x\)-x — s± (x , £i)} 
is the hyperplane of all points having the same distance between Xq and x±. The unit 
vector uj+ (xq, x\) is orthogonal to the plane H+ (xo, x\) and {x$, x\) is the oriented 
distance of that plane from the origin. 



Inversion on general domains. Suppose that ft C R™ is a bounded convex 
domain in R™ with smooth boundary dfl and denote by \Q '■ R™ ~~ ► R the characteristic 
function of 51 (taking the value one inside the domain and the value zero outside). 
Further, for any C°° function / : R™ — > R with support in fl and every xo G f2, we 
define 



(fcnf) (xo) / kn(x a ,xi) f (xi)dxi 



(1.3) 
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Fig. 1.2. Midplanc between two distinct points xo and x\ in R n . The hyperplane (xq,x\) = 
{x £ M n : a;* (xo, x\)-x = s* (xq,xi)} is the midplane betweenxi andxo, thatis, consists of all points 
having the same distance between these two points. The unit vector u>* (xq, £i) = (xi — xo)/\xi — xq\ 
is orthogonal to the plane H+ (xo, x±) and s+ (xo, xi) = (\x±\ 2 — \xo\ 2 ) / (2\x± — xo\) is the oriented 
distance of that plane from the origin. 



with 



kn (xq,xi) := < 



\xi— xo\ n r " 



2" + 1 7r"" 1 \x 1 -x r- 1 



if n is even 



if n is odd 



(1.4) 



Here and in similar situations, <9" denotes the n-fold composition of the differentiation 
operator d s - 

Remark 1.1. Since £1 is assumed to be a convex domain with C°° boundary, the 
Radon transform of xn * s a smooth function except for pairs (cj, s) G S n ~ 1 x R where 
the corresponding plane {x € R n : U) ■ x = s} is tangential to the boundary dQ. For 
two distinct points Xq,Xi inside the domain ft, the midplane H+ (xq,Xi) — {x € R" : 
(xq,x\)-x = s+ (xo,Xi)} between these points is never tangential to the boundary of 
the domain (see Fiaure UT^) and further the operators d™ andH s preserve the locations 
of singularities. As a consequence, the functions {d"H s TZxn)(^* (xq, xi) , s* (xq, x\)) 
and (d™lZxn)(u± (xq,xi) , s* (xo, xij) appearing in the nominator of the definition of 
fcn are C°° onflxil. Since additionally, \x\ — xq\ ™ +1 is weakly singular, the function 
JCnf is well defined and the operator ICci is a smoothing integral operator. 

The inversion formulas we establish in this paper are exact modulo the integral 
operator ICq. They look somewhat different in even and in odd dimensions and are 
stated in separate theorems below. 

In even dimension, our main result is as follows. 

Theorem 1.2 (Inversion in even dimension). Let n > 2 be an even natural 
number, let C R™ be a bounded convex domain with smooth boundary dfl, and let 
f: R™ — > R be C°° and vanishing outside Q. 
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Then, for every Xq G CI, 
/ fa) = (/Co/) (a: ) + 1 J ^„ "- 1 x 

OO 



V X0 ■ v x I ± ■ ^ drdS (x) , (1.5) 

Jan Jo r z — \x — x \ 

/ (ar ) = (/Co/) (ao) + ^ X 

■ (x - / — : k drdS W ■ I 1 - 6 ) 

an Jo r z — \x — sco | 

In feoi/i formulas, the inner integration is taken in the principal value sense, dS is 
the usual surface measure, v x is the outward pointing unit normal to CI, and /Co 
is the integral operator defined by il.3\) . Moreover, D r := (2r)^ 1 d r denotes 

differentiation with respect to r 2 and V Ko • the divergence with respect to xo ■ 
Proof. See Section [2] □ 

In odd dimension we have the following corresponding result. 

Theorem 1 . 3 (Inversion in odd dimension) . Let n > 3 be an odd natural number, 
let CI C K" be a bounded convex domain with smooth boundary dCl, and let f: TSL n — > K 
be C°° and vanishing outside of CI. 

Then, for every Xq G CI, 



f (xo) = (/Co/) (xo) + K 4—1 — x 



(-l)^- 3 )/ 2 ^ 
4tt t1 



V X0 -/ Vx (V r r 2 r n - 2 Mf)(x,\x-x \)dS(x) , (1.7) 
Jan 



f (xo) = (ICnf) (xo) + V ' , n _ x x 



( _l)(n-3)/2 a;n _ i 
47T" 

v x -(x -x)(d r V?- 2 r n - 2 Mf)(x,\x-x \)dS{x) . (1.8) 

an 

Here /Co, v x , V l0 , dS, and T> r are as in Theorem M.'SX 
Proof. See Section HI □ 

Both, Theorem 11.21 and Theorem 11.31 will follow from corresponding statements 
for the inversion of the wave equation in even and odd dimensions, which we shall 
establish in the following sections (see Theorems 13.11 and I4.1j) . 

Exact reconstruction for elliptical domains. In the case that CI is an ellip- 
tical domain, we show that the integral operator /Co vanishes exactly and therefore 
Theorems 1 1 . 21 and 1 1 . 31 provide exact reconstruction formulas for ellipsoids. After trans- 
lation and rotation, we may assume that the elliptical domain takes the standard form 

CI := G R n : \A^x\ 2 < l} , 

where A G W ixn is a diagonal matrix with positive (possibly distinct) diagonal entries. 
Obviously, a ball is a special case of an elliptical domain where all diagonal entries of 
A coincide and are equal to the radius of the ball. 
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For elliptical domains we have the following exact inversion formulas. 

Theorem 1.4 (Exact inversion on elliptical domains). Let 57 C R ra be an elliptical 
domain. Then /Co/ vanishes identically on CI. In particular, for any smooth function 
f : 1" — > R with support in CI and every xo G CI, the following hold: 

(a) If n is even, then 

/(*») = ^ X 



v / f°° (rV n r - 2 r n - 2 Mf){x,r) a 

Vx • / v x \ ± ■ ^ drdS(x). (1.9) 

on Jo r z — \x — x \ 



(b) If n is odd, then 

/(*») = ^ x 



V X0 - v x (V?- 2 r n - 2 Mf)(x,\x-x \)dS(x). (1.10) 
Jan 



Here v x , V l0 , d5* ; and T> r are as in Theorem 
Proof. See Section [5] □ 

Since Knf = 0, Equations (|1.6j) and (|1.8|) provide alternative exact inversion for- 
mulas for elliptical domains. We omit formulating separate statements. Further, by 
taking limits one can easily establish exact inversion formulas like (|1.9p and p,10[) for 
certain unbounded domains, such as for elliptical cylinders. We again omit formulat- 
ing such generalizations. 

1.2. Relations to previous work. Exact back-projection type inversion for- 
mulas for recovering a function from spherical means with centers on the boundary of 
a ball have been discovered quite recently in [7J [HI HH HS1 [23] ■ In [23] a formula has 
been found for n — 3, which has later been generalized to arbitrary dimensions in |12j . 
In odd dimension these formulas coincides with our Equation (jl.101) (which, however, 
holds for the more general case of elliptical center sets) . A different set of exact re- 
construction formulas has been derived in [8] for odd dimensions and in [Jj for even 
dimensions. In [9] relations between the different formulas have been investigated for 
dimensions n — 2 and n = 3. None of these papers considers the case of more general 
domains. In [13] reconstruction formulas of the back-projection type have been found 
for certain polygons and polyhedra in two and three spatial dimensions. 

Formulas that recover a function from spherical means with centers on the bound- 
ary of an elliptical domain in arbitrary dimension have been obtained in [18] . The 
derived identities as well as the method of proof are different from ours. Our results 
are, however, closely related to ones of pUHH]. Actually, the present article general- 
izes the result obtained for n = 2 in [10] and for n = 3 in [M] to the case of arbitrary 
spatial dimension. 

1.3. Outline. The main aim of the following sections is the proof of Theo- 
rems II .21 11.31 and 11.41 To that end, we first derive an auxiliary identity for the wave 
equation in Section [2] (see Theorem I2.2[) . Subsequently, in Section [3] we shall prove 
Theorem 11.21 and in Section [4] we establish Theorem 11.31 In these sections, we also 
derive corresponding statement for recovering the initial data of the wave equation 
from the solution on the boundary of an arbitrarily shaped domain. These results, 
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which are also of interest in its own, will be presented in Theorems 13.11 and 14.11 be- 
low. In Section [5] we consider the case of elliptical domains, where we show that the 
operator ICq vanishes identically and therefore we establish the exact reconstruction 
formulas stated in Theorem 11.41 The paper concludes with a discussion in Section [6l 

2. Auxiliary results for the wave equation. Let ft C R™ be a bounded 
convex domain in R™ with smooth boundary dtt and let /: R™ — > R be a smooth 
function with support in tt. Consider the following initial value problem for the wave 
equation 

' (<9 t 2 - A x ) p(x,t)=0, for (x, t) e R" x (0, oo) 

p(x,0)=f(x) , forxeR" . (2.1) 

d t p (x, 0) = , for x e R" 

Here dt denotes differentiation with respect to the temporal variable t £ (0, oo) and 
A x is the Laplacian in the spatial variable x € R™. To indicate the dependance of the 
initial data we will also write the solution of (|2.1|) as p = Wf. 

According to the well known explicit formulas for the solution of (|2.1j) in terms of 
spherical means, recovering a function from spherical means is essentially equivalent 
to the problem of recovering the initial data in (|2.ip from values of the solution on 
<90 x (0, oo). In this section we derive a basic result for the wave inversion which 
in the following sections will be applied to derive the results for the inversion form 
spherical means presented in the introduction. 

2.1. Outgoing fundamental solution. Throughout the following we denote 
by G (x, t) the outgoing fundamental solution (or free space Green's function) of the 
wave equation, that vanishes on {t < 0} and satisfies the equation 

(<9 t 2 - A x ) G (x, t) = S n (x) 5 (t) , for all (x, ijeff'xl. 

Here and in the following, S n and S denote the n-dimensional and one-dimensional 
delta distribution, respectively. 

Remark 2.1. Note that also in the even- dimensional case we consider G (x , t) as 
a distribution on R™ x R and not as a classical junction. The argument (x, t) in G (x, t) 
is only a formal notation ( and, in particular, does not mean the point-wise evaluation 
of G at (x, t) ) that indicates the existence of a spatial and a temporal variable where 
these distributions act on. Derivatives of the fundamental solution, like dtG(x,t), 
will always denote distributional derivatives. 

With the outgoing fundamental solution of the wave equation, the solution of the 
initial value problem (|2.1j) can be written as 

p(x,t)= I (dtG(x-Xut))f(xi)dxi, for all (x, t) e R" x (0, oo) . (2.2) 

Note that throughout this paper integrals like the one on the right hand side in (|2.2I) 
will always be understood in the week sense. Hence, the identity (12. 2\i actually means, 
that we have 

f [ v(x,t)p(x,t)dxdt= [ (d t G Xl ,v)f(x 1 )dx 1 (2.3) 

for any smooth test function v : R™ x R — >• R. Here ( ■ , ■ ) denotes the duality pairing 
between a distribution and a test function, p (x, t) is extended with zero on {t < 0} 
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and G X1 (x,t) = G (x — x±,t) denotes a shifted copy of the fundamental solution G 
translated by X\ € W 1 . 

After inserting the known explicit expressions for the outgoing fundamental solu- 
tion G, Equation (|2.2[) can be rewritten in terms of spherical means of the function /. 
The explicit solution formulas differ in even and in odd dimensions and will be stated 
in the following Sections [3J and HI where we study the even and the odd dimensional 
case separately and in more detail. 

2.2. Universal backprojection. For any smooth function v. dft x (0, oo) — > K 

we define 

(Bov)(xo) :=2V X0 • / v x [ G (x - x , t) v (x,t) dtdS (x) , for^Gfi. (2.4) 
Jan Jr 

Note that in the even dimensional case, the function v (x, t) needs some decay as 
t — > oo in order that the integral (Bqv) (xq) is well defined. This is certainly the 
case if v is the restriction of the solution of the initial value problem (|2.ip for some 
smooth initial data / supported in f2, which happens in all instances where we apply 
the operator Bq. In three spatial dimensions (and in a slightly different form), the 
inversion integral (|2.4p has been introduced to photoacoustic tomography in [23 . 
According to the notion of [23j . we call Bq the universal backprojection operator. 

In [23] it has been shown that the identity B^Wf = f holds for the case that 
£1 is ball in three spatial dimensions and that / vanishes outside f2. This exact 
reconstruction property does not hold for general domains. However, for arbitrarily 
shaped domains in arbitrary dimension we have the following result: 

Theorem 2.2. Let fl C M n be a bounded convex domain with smooth boundary. 
Then, for any C°° function f : M" — > K with support in Q and any xo € ^2, we have 

f (x ) = (B n Wf) (x ) + f f ( Xl ) *4 2) (x , x 1 )dx 1 , (2.5) 
Jn 

with the distributional kernel 

k^ ) (x ,x 1 ):=(V X0 +V Xl ) 2 ( f (d t G(x~ Xll t))G(x-x ,t) dtdx. (2.6) 

Jn Jr 

Here (V Xo + V Xl ) 2 is s shorthand notation for the operator (Vx + V 2 , 1 )-(V a;o + V Xl ) 
and the equations (2. 5)) and \2. 6]) have to be read in the weak sense. 

Proof. According the the well known Kirchhoff integral representation, the solu- 
tion p = W f of the initial value problem (|2.1[) satisfies 

f{x )= / Vx - \ {V x p(x,t))G(x- x ,t)dtdS(x) 
Jan Jr 

I v x - I p(x,t) (\J X G (x - x ,t)) dtdS (x) . 
Jan Jr 

Now inserting the identity GV x p — V x (pG) — pV ' X G in the first term followed by 
an application of the divergence theorem, and using the relation V X G (x — XQ,t) = 
~V Xo G (x - x ,*) yield 

f(xo)= [ [ A x {p{x,t)G(x-x ,t))dtdx 
Jn Jr 

+ 2 V X0 • / u x I p (x, t) G (x - x , t) dtdS (x) . (2.7) 
Jan Jr 
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According to the definition of the second term equals (SfiW/) (xo)- After insert- 
ing the representation (|2.2[) for the solution p = W/ of the initial value problem (|2.1[) 
and applying the relation X7 X G (x — Xi,t) = —\7 X .G (x — Xi,t) with i = 0, 1, the first 
term in (|2.7|) is seen to take the form L fc s ^(x , x\)f {x\) dxi, with fc^ defined by 
Equation (|2.6I) . This concludes the proof of Theorem 12.21 □ 

Remark 2.3. In the proof of Theorem \ as well as in the following deriva- 
tions we formally operate with distributions as they were classical functions. These 
computations can be made more rigours by writing down all equalities in the weak 
sense (similar as done in UOf for the two dimensional case) and then using classical 
integral calculus. However, the calculus using distributions used in the present paper 
seems to be more intuitive and easier to follow and has also been used in JH 1 1$ for 
the derivation of inversion formulas in three spatial dimensions. 

3. Inversion in even dimension. Now let n > 2 denote an even natural num- 
ber. In this section we derive explicit formulas for the wave inversion in even dimension 
and then apply these results for establishing Theorem 11.21 In even dimension, the 
outgoing fundamental solution of the wave equation takes the following explicit form 

.p(„-a)/2 x{*M*l 2 >o} on{t>0} 



G(x,t) = { 2W2 * , (3.1) 

[0 on {t < 0} 

where V t = (2t) 1 dt denotes differentiation with respect to t 2 , and x {t 2 ~ \ x \ 2 > 0} 
is the characteristic function of the set all points (x, t) £ R™ x R with t 2 — \x\ 2 > 0. We 
emphasize again that in Q3.1[) and in similar situations all derivatives are understood 
as distributional derivatives. 

We now have the following result for recovering the initial data of the initial value 
problem (|2.1I) from the restriction of its solution to dtt x (0, oo). 

Theorem 3.1 (Wave inversion in even dimension). Let n > 2 be an even natural 
number, let £1 C M n be a bounded convex domain with smooth boundary, and let 
f : M. n — > R be a C°° function that vanishes outside of CI. Then, for every xq G Q, 



f(x ) = (ICnf) {x ) 



(-1) 



(n-2)/2 



7r«/2 

Vx a ■ v x ~ , = dS (x) , (3.2) 

Jon J\x-x \ Jtft-\x-x \ 2 

/ uM)/2 
f{x Q ) = {Knf) {xo)+ K -^j- 2 x 

v x ■ (Xq -x) = dS (x) . (3.3) 

•^-^l \jt 2 - \x ~ X Q \ 2 

Here ICq, v x , \7 xo , and dS are as in Theorem ] 1. 6 A 

We proceed this section by first deriving Theorem 13.11 and then establishing 
the corresponding result for the inversion from spherical means in even dimensions 
(namely Theorem 11.21) . 

3.1. Proof of Theorem 13.11 According to Theorem 12.21 we have to show that 
the kernel fc^ 2 ' defined in (|2.6|) is equal to the kernel fcn defined in (|1.4p . and that 
(BnWf) (xq) can be written as the integral term in (|3.2I) as well as the one in 
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Let us start by showing that = fen, that is, 
(V ao + V.J 2 f [ (d t G(x-x u t))G(x-x ,t)dtdx 

I I 2 I I 2 

where = (ar ,£i) = \*]ZZ° \ and s * = s * ( x 0,%i) = 2 | li-so | are as in C2])- 

To show (|3.4| . for any two given points xo ^ x\ G f2, we write i?o := \x — xq\ and 
i?i := |as — affx | - Then, using the explicit expression ()3.1| for the fundamental solution 
of the wave equation in even dimensions, the inner integral on the left hand side of 
(|3.4p evaluates to 

(d t G (x — x\,t)) G (x — xo,t) dt 

= J_ r ( d x{t 2 -Rj>o} \ (w - a)/a x {t 2 - fl§ > o} 



— r ^ P (n-2)/2 X ~ ^ > 0} ^ (n-2)/2 X {t 2 - fl§ > 0} 
J_ c (n-2)/2^(„-2)/2 f X{t 2 -i? 2 >0} \ X {f 2 -.Rg>0} 

4^" fii % 7o v * v^ 3 ^! / 

i-T 



(„-2)/2 („-2)/2 j. /• 

' Kl % r^cx, Hl J m: 



tdt 



2 



For T > max{i?0j -Ri}> the above integral on the right hand side computes to 
rT tdt 



1 



= ^^-^ + ^-^--1^(1^-^!) . (3.5) 

After applying the operator and letting T — > oo, the first term vanishes. Now 
let $ (s) = 1/s denote the principal value distribution <p H> lim^o Jr\(_ £ e ) <^(s)s _1 ds. 
Recalling the definitions of Ro and R\ then implies 



| (ftG (x- Xl ,t))G(x-xo,t)dt = -^ ■D%- 2)/2 V < £ - 2 ' " 



1 



( i W 2 

= ^T" $( "" 2) (l* - ^o| 2 - |* " xx| 2 ) . (3.6) 

Here and in the following denotes the z/-th distributional derivative of $ for some 
integer number v > 0. 

I 2 I 1 2 

For the following recall that oj* = 11-3:0 and = ~ , and write any point 
x G K™ in the form x = .sw t + y with s6K and yJ-W. We then can compute 

\x - x \ 2 - \x - a;i| 2 = (s - cj* • x ) 2 - (s - cj* • xi) 2 

= 2uj+ ■ (x - xi) (s* - s) = 2 |xi - aco| ( s - s *) • (3-7) 
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Together with Equation (13.61) and the definition of the Radon transform 1Z, this further 
implies 



(d t G (x — x\,t)) G (x — xo, t) dtdx 

= —T^r- / / Xn (aw* + y) $ ( "- 2) (2 |n -x \{s- a*)) dy ds 
= { -^- J r $( "~ 2) ( 2 l*i " *o| (s ~ **)) ^ ^ Xo + 2/) dy) ds 

(_-t\n/2 

= M^" / $ ( "" 2) ( 2 l^i " *o| (s ~ «*)) (ftxn) K, a) ds . 



Now using the chain rule, integrating by parts n — 2 times, recalling the definition 
of the principal value distribution $ (s) = 1/s, and noting that the Hilbert transform 
H s is defined as the convolution with 7r _1 $ imply 



(dtG (x — xi,t)) G (x — xo, t) dtdx 
(-l) n/2 



47r«2™- 2 |n - a; |" _2 



(G>r 2 $ (2 \xi - x \ (s - s*))) (ftXn) (w*, a) ds 



f-lV* /a 

' (9r 2 ^Xn)(^,s)d S 



7r n 2™ |a;i - x | n 2 Vr 2 |n -x \(s- s*) 
(-1) ( "^ 2)/2 /" 1 



2 n+l 7r n _ ^j" 1 J R Si< - s 



(df- 2 n X n)(^,s)ds 



f_i\M)/2 



2™+ 1 7T"- 1 n - 3f0 



It remains to apply the operator (V Xo + V^) to the last expression. By using 
the relation (Va, + V Kl ) s* = (n — xq) / |n — xq\, the result is 



(V* +V Xl ) 2 / / (atG?(aj-a;i,t))G(a;-a;o,*)dtdx 



f_i\M)/2 



2 n+l 7r n-l ^ _ 



This is the equality claimed in (|3.4[) . 

Now recall the definition of Bq (see Equation (|2.4[l ^l as well as the explicit rep- 
resentation (|3.1[) for the fundamental solution of the wave equation in even dimen- 
sion. Further notice that for any integer the formal L 2 adjoint of T>" is given by 
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(2?j) = (— 1) W"t . We therefore can compute 



(BaWf) (s ) 

1 / [°°„( n -2)/2 ( X{t 2 > \X-X \ 2 } 



(-i) ( "' 2)/2 v . f „ r *{* 2 >i"oi 2 } 



r™/2 



x 



an J ° Ji 2 - |a;- sol 



\Jt 2 - \x- XoT 

^ p („-2)/2^i w ^ ^ dtdS*(a;) 



T n/2 



V X0 - / v x / v 1 . -LL± ' ; dtdSfr) 



F - x a | 



2 



In view of Theorem 12.21 and Equation (|3.4p . the above expression for (BqW/) (xo) 
yields the first identity in Theorem 13. 11 formula (|3.2j) . 

Finally, we verify the second identity in Theorem 13.11 Interchanging the order 
of differentiation in the last displayed expression for (BnWf) (xo) followed by one 
integration by parts yields 

( _ 1) ( n -a)/2 7r -n/a {BnWf) [xq) 



on 



/',■(.«■„-.<■) / ih | X{t 2 > \X X \ 2 } | V {n-2)/2 t - lw/{r ndulS{r) 

t 2 - \x - x \ 2 



90 ^|x -x| ^ t 2 _ \ x _ ^J 2 

This shows equation (|3.3I) and concludes the proof of Theorem 13.11 

3.2. Proof of Theorem 11.21 Recall the formula (|2.2[) for the solution of the 
wave equation (|2.1[) as well as the explicit expression (|3.1[) for the fundamental solution 
of the wave equation in even dimensions. After introducing polar coordinates around 
the center x € dfl we can write 

(yVf)(x,t)= [ {d t G{x-x u t))f(x 1 )Ax 1 



n 



1 f I Wn-2)/2X{* 2 -F-3=li 2 >0} 



2tt"/ 2 



y i 2 - \x - x\Y 



' r" L Mf{x,r) \d t V\ " -4=4 dr. 



2W 2 Jo I 



After multiplying the last displayed equation with G {x — xq, t), integrating over 
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the time variable and using the shorthand notation i? = \x — Xq\ we obtain 



/ G(x-x ,t) Wf (x,t) dt 
Jr 



u n -i f°° ( p(tt-2)/a X {t 2 > ^0} 



4tt" 



r^M/^r) L['^)/^{ t2 > ^ drdt 
o V v - r 2 ) 



T n - X Mj[x,T) X 



(I 



" 4tt" 



r"- 1 ^ / (z, r) I3 ^- 2 )/ 2 I3 (™-2)/2 x 



T 2idf 



D r lim / , dr . 

\ T->oo J max{Ro , r} y/j* - i?2^2— 72 J 

The inner integral has already been computed (see Equation (|3. 51) ) and shows 
X>r I™ / , = 5 -j = (r 2 - i? 2 ) , 



max{ii ,r} ^ 2 - i^V^T 2 " ^ 2 - R t 







with $ denoting the principal value distribution P.V. \/s. After recalling that 2? r 
denotes differentiation with respect to r 2 and that the formal L 2 adjoint ol V" is 
given by (V")* — {— 1)" rV v r r~ x we obtain 

G(x-x ,t) Wf (x,t) dt 

, ,x(n-2)/2 -oo 
( — J- J W n _l / 



y r n_1 M/ (sc, r) I?; 1 " 2 * (r 2 - |a: - x | 2 ) dr 



47T n 

(_ 1 )("-a)/2 Wn _ i ,00 (r2?»-2r n - 2 >t/)(x,r) 



4tt" Jo r 2 - |x - jc | 



2 



dr . 



Finally, by using the definition of Bn (see Equation (|2.4p ) and inserting the identity 
just established we obtain 

(BnW/) (a*) = 2 V X0 • / v x [ G{x- x , t) Wf (x, t) dtdS (x) 
Jan Jr 



z " -/an Jo r z — \x — xo\ 



According to Theorem 12.21 and Equation (|3.4[) this shows the first inversion for- 
mula (|1.5[) in Theorem 11.21 

It remains to verify Equation (ll.6[) . This equation, however, is an easy con- 
sequence of the identity (|1.5[) just established. In fact, interchanging the order of 
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integration and differentiation in (|3.8[) and integrating by parts yields 

( _,,(n-2)/2 

(BnW/) (x ) = ( 1} - x 

i/ a • (ac ~ a;) / : k drdS \ x ) • 

ao Jo r z — \x — xq\ 

Again, according to Theorem 12.21 and Equation (|3.4p the last displayed equation im- 
plies (|1.6p and concludes the proof of Theorem 11.21 

4. Inversion in odd dimension. Now let n > 3 be an odd natural number. In 
this case, the outgoing fundamental solution of the wave equation is given 



p(n-3)/a 



8 (i 2 - |a;| 2 ) on {t > 0} 



G(x,t) = < 2^1771 t - V ,-, y — i > (41) 

(0 on {t < 0} 

Here, as usual, the operators V t = (2i) _1 d t denotes the distributional derivative with 
respect to the variable t 2 . 

We have the following counterpart of Theorem 13. II for odd dimensions. 

Theorem 4.1 (Wave inversion in odd dimension). Let n > 3 be an odd natural 
number, let CI C M n be a bounded convex domain with smooth boundary, and let 
f : 1™ — > R be a C°° junction that vanishes outside of CI. 

Then, for every Xq £ CI, we have 

f (xo) = (fcnf) (xo) 

+ 2n (n-i)/2 V *° ' j m v * (^ ( "" 3)/2 ^ 1 W/) (x, \x - x Q \) dS (x) , (4.2) 

/ (X Q ) = (/Co/) (x ) 

+ 27r (n-i)/2 j m v * ■ ^ ~ x ) (dM n - 3)/2 t^Wf) (x, \x-x \)dS (x) . (4.3) 

Here, again, /Co, v x , V xo , and dS are as in Theorem ] l.'A 

We proceed with this section by first establishing Theorem 14.11 and then deriving 
the formulas in Theorem IL3J as corollaries of it. 



4.1. Proof of Theorem 14.11 Similar to the even dimensional case we apply 



Theorem 12.21 and verify that the kernel k^ ' defined in (I2.6[) is equal to the kernel fco 
defined in (|1.4j) . and that (BnWf) (xq) can be written as any of the integral terms in 
Equations P~2"j) and P~3l . 

We first show that fc = ka, that is, 

(V Xo +V ai ) 2 / / (d t G(x-xi,t))G{x-x ,t)dtdx 



(_-\\(.n-l)/2 

r (^ Xn )(w*,**) , (4.4) 



2 n+l 7r n-l \ Xl _ ^ 

where u;* = (x ,a;i) = [f^Efg] and s* = s* (x ,a;i) = ^73^] are as in ((Ol)- 
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With the notation i? := \x — x \ and i?i := |a; — n|i the representation (|4.1I) for 
the outgoing fundamental solution of the wave equation in odd dimensions yields 

(d t G (x — n,i)) G (x — xo,t) dt 

i r°° / \ 

4 J^ I (trtr 1)/2 5 f r\)) rtr* )/2 5 f - Ri) ^ 

2 - £ (tV%-^ 2 6 (t 2 - Rl)) V^S (t 2 - i? 2 ) dt 



= <- 1)/2 <- 3)/2 Jf S (t 2 - fl?) 5 (t 2 - i? 2 ) 2tdi 



/ 1 x(n-3)/2 



4^" 



— 5 ( "- 2) (ix-xol 2 -^-^! 2 ) . 



Here and in the following denotes the ^-th derivative of the one-dimensional delta 
distribution for some integer number v > 0. 

Now recall the definitions = Xl ~ x ° and s* = Hrr ~^"| and write a; = sw t + w 

|X1— Kol 2|£Ci-K | w 

with s£l and y-Lw. We then have (see Equation (13. 7[1 ) 

|x — xo| 2 — |a: — a;i| 2 = 2 |n — aco | (s — s*) . 

This implies 

(dtG (x — Xi,t)) G (x — xq, t) dtdx 

X n (slj* + y) <5 ( "~ 2) (2 |n - x \ (s - s*)) dyds 
,5("- 2 ) (2 |n - x | (s - s*)) (Y Xn (sw* + y) dy) ds 
<5("" 2 ) (2 |n - i | (« ~ «*)) (^Xn) K, *) ds . 



(- 


.!)(«- 


-3)/2 




4tt™- 


-1 


(- 




-3)/2 




47T™- 


-1 


(- 


.!)(«- 


-3)/2 




4tt™- 


-1 



Integrating n — 2 times by parts yields 

(dtG (x — xi,t)) G (x — xo, t) dtdx 



47T"- 1 2«- 2 |ari-a | 



n-2 



J (2 In - i | (« - **)) (d^Uxa) (w*, s) ds 



C_ n («-l)/2 

— (9»- 2 ^)(^,s.) . 



2 «+l 7r n-l | Xl _ Xo |' 

As in the even dimension case, after application of (V Xo + V Xl ) 2 this yields (|4.4j) . 
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Next, note that the fundamental solution (|4.1[) in odd dimensions may be rewrit- 
ten in the form 1/ {A / n^ n ^ 1 ^ 2 )T>[ n ^ 3 '^' 2 t~ 1 8 (t — \x\). Consequently, by the definition 
of Bn (see Equation (|2.4[1 ). we have 

(B n Wf) (x ) = 2 V X0 ■ [ v x f G(x- x , t) Wf {x, t) dtdS (x) 
Jon J& 

= 2 ^(n-l)/2 V *o ' J V * (^ (n_3)/2 *- 1 W/) (X, \X - X Q \)dS (X) . 

In view of Theorem 12.21 and due to Equation (|4.4[) this implies formula (|4.2j) claimed 
in Theorem 14. II Finally, carrying out the differentiation under the integral yields 

(B Q Wf) (s ) 

= 27r (n-l)/2 j m v * ■ ^ ~ X ) (M n ~ 3)/2 * -1 W/) (*> I* " *o|) (z) . 
This shows that also identity (|4. 3[) holds and concludes the proof of Theorem 14.11 

4.2. Proof of Theorem ll.3l Inserting the expression (|4.1j) for the fundamental 
solution of the wave equation in odd dimensions (see Equation (|2.2[1 ) shows that the 
solution of the wave equation (12 . If) can be written as 

= 47r 7"ri)/2 ^ n - 3)/a r n - a ^/(i,r) . 
Together with the definition of 2?n (see Equation (|2.4t ) this yields 
(BnW/) (x ) 

= 2^(«-D/2 ' ^ n (^ n - 3)/2 »- 1 W/) (x, |s - a;o|) dS (s) 



^V I0 ./ ^ (l?("- 3 )/ 2 r- 1 ^"- 3 )/V"- 2 X/) (x,|x-x |)d5(x) 
^ Jan v 7 



8 

'an 



4 _„_ a V xo - /_ _v x {V?- 2 r n - 2 Mf)(x,\x-x \)dS(x) 



According to Theorem UM and Equation f|4.4[) this yields ()1.7[) . 

It remans to establish the second formula in Theorem ll.31 namely Equation (|1.8p . 
To that end, one simply carries out the differentiation in the last displayed formula 
for (BnWf) (so), which yields 

(BnWf) (s ) = V X0 • / i/ B (I?r 2 r"- 2 M/) (s, |s - x \)dS (x) 

47r Jan 



v x ■ (s - s) (d r V?- 2 r n - 2 Mf) (x, \x - x \) dS (x) 

an 



47T™- 1 

This however yields formula (|1.8j) . 

5. Exact inversion for elliptical domains. Let A — diag (ai , . . . , a n ) be a 

diagonal matrix in R" XTl with entries ctj > and consider the elliptical domain 



n 



<js G R n : \A~ l x\ 2 < l} 



In order to establish the exact inversion formulas of Theorem 11.41 it is sufficient to 
show that JCnf = 0. This will be done by first verifying that the kernel vanishes for 
the special case that the domain is a ball and then applying a linear transformation 
to the ball to establish the result for general case. 
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5.1. Special case: spherical domains. Let B :— {x £ M." : \x\ < 1} denote 
the unit ball in R™ centered at the origin. Then, elementary geometry shows that the 
Radon transform of xb is given by 

n XB (w,s) = < v > 1 1 . , (5.1) 

I otherwise 
where = ""~ 2 denotes the volume of the unit ball in R™ _1 . 

n 1 n— 1 

Odd dimension. If n is odd, then (15.11) shows that 7Z\b (w, s) is a polynomial of 
degree n— 1 on {|s| < 1}. Since |s* (xo, £i)| < 1 for any two distinct points x$, x\ € B, 
this yields 

( _ 1)( n-l)/2 (0nft ) (a . 0) (a . 0) a;i)) 

Fi - x | 

This implies that we also have /Cs/ = and, according to Theorem ll.31 this shows the 
inversion formula (|1.10|) stated in Theorem [L4] for the special case that the considered 
domain is a ball in odd dimension. 

Even dimension. If n > 2 is an even natural number, then the identity JCsf = 
is slightly less obvious. In this case, we first note the following identity (see, for 
example, [201 Table 7.3, Number 13]) 

(H s (sct>)) (s) = s {n s 4>) (§)-[<(> (*) ds (5.2) 

JR 

satisfied by the Hilbert transform and some function <fi: M. — > R. Further, Equa- 
tion (|5.ip shows that we have the relation (TZxb) ( w , s ) = Pn-2 ( s )0i/2 ( s )> where 
Pn-2 ( s ) is a polynomial of degree n — 2 and 4>\/2 (s) '■= \/max{0, 1 — s 2 }. Applying 
the identity (I5.2[) repeatedly, thus yields 

{H s TZxb) (w, s) = Q„_2 (s) (Hs^l/2) (s) + Qn-3 (s) , 

for certain polynomials Q n -2 and Q n -3 of degree n — 2 and n — 3, respectively. The 
Hilbert transform of (f>i/ 2 (s) = ^/max{0, 1 — s 2 } is known and given by (see, for 
example, gj Table 13.11]) 

(^1/2) 0) = -s + sign(s)x{|s| > 1} Vs 2 - 1 for all seK. 

In particular, (% s (j>i/2) (s) is a linear function on {\s\ < 1} and therefore the product 
Qn-2 ( s ) ( s ) is a polynomial of degree n — 1 on {|s| < 1}. Noting again that 

|s* (xo,a;i)| < 1 for any two distinct points xq,x± € B, we therefore conclude 

, / \ (- 1 ) ( "~ 2)/2 {d™H s K XB ) (u* (x ,xi) , (x , xi)) 

fcs (aso.ari) = 9n+ i_ n -i ; = for ^ ^ x e B . 

Fi - sol 

This implies that the identity /C_b/ = also holds in even dimension. In view of 
Theorem 11.21 this establishes inversion formula (|1.9|) in Theorem 11.41 for the special 
case that the considered domain is a ball in even dimension. 
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5.2. General case: elliptical domains. Now let 57 = {x £ W l : |A _1 2;| < l} 
be an elliptical domain where A = diag (ai, . . . , a n ) is a diagonal matrix with positive 
entries that possibly differ from each other. We then obviously have the identity 
Xn {%) = Xb (A~ 1 x^, where B C E™ is the unit ball considered in the previous 
subsection. Therefore, the known relation between the Radon transform of a function 
ip and the Radon transform of the function x i— > ip implies that 

K X n («, s) = K X B ( , for all ( w , a) e x R , (5.3) 

From (15.31) we conclude that 

(d?K X n) (w, a) = (^ X b) f TTT- TTt) if " is odd . 

\Alj\ + \\ A ^\ \ Auj \J 

(Wxn) (w, ») = (WW) f 7X7- TXl) if « is even ■ 

According to the special case considered in Subsection l5 . 1 I this shows that kn (xq, x\) = 
for all xo ^ x\ £ £1 and hence that Knf — 0. In view of Theorems 11.21 and 11.31 
this establishes the exact reconstruction formulas (jl.9|) and (| 1 . 10[) for the inversion of 
spherical means on elliptical domains in arbitrary spatial dimension. 

6. Discussion. Many medical imaging and remote sensing applications aim for 
recovering a function from spherical means centered on a set of admissible receiver or 
detector locations. In the case that the center set is an infinite hyperplane explicit 
formulas of the back-projection type for recovering a function from spherical means 
are known since the mid 80s (see [HE]). In the case that center set is a spherical 
or cylindrical surface such type of formulas have been derived about 20 years later 
in HI [HI [23] . All these geometries have rotational and/or translational invariance 
and seem well adapted to the inversion from spherical means. It therefore has been 
believed by many researchers that such exact back-projection type inversion formulas 
may only exist for those invariant geometries. 

Very recently, in |1Q[ 114] explicit exact inversion formulas of the back-projection 
type for inverting the spherical mean transform with elliptical center sets in dimen- 
sions n = 2 and n = 3 have been derived (see [18] for a different formula for ellipsoids 
in arbitrary dimension and 13{ for reconstruction formulas for certain polygons and 
polyhedra). Moreover, in [TU1 E] it has been shown that the same formulas may be 
applied when the center set equals the boundary of an arbitrarily shaped smooth con- 
vex domain 57, in which case these formulas recover the unknown function modulo an 
explicitly computed smoothing integral operator ICq,. In the present paper we gener- 
alize these results to the case of arbitrary spatial dimension. We have further shown, 
that the operator /Co vanishes for elliptical domains which yields exact inversion for- 
mulas in these cases. However, as can be readily verified by using Equation (|1.4I) . the 
operator JCq does not vanish for general domains. Actually, these results give an affir- 
mative negative answer to the question whether the universal back-projection formula 
of Xu and Wang [23] (introduced there for the case n = 3) is exact for general do- 
mains. This negative result, however, does not imply that a different back-projection 
type formula may provide exact reconstruction for general domains. Nevertheless, it 
is believed by the author that such a truly universal reconstruction formula for the 
spherical mean transform does not exist. 
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